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We calculate quarkonium spectral functions in a quark-gluon plasma using a potential model based 
on full QCD lattice calculations of the free energy of static quark-antiquark pair. We estimate the 
binding energy and the thermal width of different quarkonium states. The estimated upper limit for 
the dissociation temperatures is considerably lower than the ones suggested in the recent literature. 



PACS numbers: ll.15.Ha, 12.38.Aw 

One of the most important features of the quark gluon 
plasma is the screening of static chromo-electric fields 
It has been argued that screening above the decon- 
finement temperature Tc is strong enough to lead to the 
dissolution of the J/ip state, which can then signal quark- 
gluon plasma formation in heavy ion collisions The 
fact that the J/V'-suppression pattern observed in SPS 
and RHIC experiments is still not understood 0], serves 
as a motivation for this work. In particular, our aim is 
to relate quarkonium suppression to deconfinement and 
color screening. 

Because of the large quark mass m = rricb 3> Aqcd, 
the velocity v of heavy quarks in the bound state is small, 
and the binding effects in quarkonia at zero temperature 
can be understood in terms of non-relativistic potential 
models [3]. More recently, the potential has been de- 
rived from QCD using a sequence of effective field theo- 
ries: Non-relativistic QCD (NRQCD), an effective theory 
where all modes above the scale m are integrated out, and 
potential NRQCD, an effective theory in which all modes 
above the scale mv are integrated out Q. 

Based on the success of the potential model at zero 
temperature, and on the idea that color screening im- 
plies modification of the inter-quark forces, potential 
models have been used to try to understand quarko- 
nium properties at finite T To discuss properties 
and dissolution of quarkonium states at finite temper- 
ature spectral functions have to be considered. Using 
lattice QCD, charmonium correlators have been calcu- 
lated, and the corresponding spectral functions have been 
extracted u sing the Maximum Entropy Method (MEM) 
0, B B E EJ- The MEM at zero temperature can re- 
construct the basic features of the spectral functions: the 
ground state, the excited states, and the continuum 11 1. 
At finite temperature, however, the extraction of the 
spectral functions becomes difficult, because the length of 
the Euclidean time is limited by the inverse temperature. 



In a recent paper [I4I we calculated the spectral func- 
tion using a potential model for the non-relativistic 
Green's function. This detailed analysis shows that spec- 
tral functions calculated in a potential model combined 



with perturbative QCD can describe the available lat- 
tice data on quarkonium correlators at zero and at finite 
temperature in QCD with no light quarks [l2|- Char- 
monium states were found to dissolve at temperatures 
significantly lower than quoted in lattice QCD studies. 
In the present work we extend our analysis to the case of 
real QCD with one strange quark and two light quarks. 
Furthermore, using lattice results on color screening we 
derive upper bounds for the dissociation temperatures 
of various quarkonium states in the quark-gluon plasma 
which are considerably lower than previous estimates. 

In the energy domain of the resonances and continuum 
threshold the spectral functions are calculated from the 
non-relativistic Green's functions as 
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aiu) = ^i^ImV • V'G"'-(f, ?, E)\,^~,^^^ , (2) 

for S'-wave, and P-wave quarkonia, respectively. Here 
E = u! — 2m , and the number of colors is iVc = 3. The 
non-relativistic Green's function satisfies the Schrodinger 
equation + V{r) - e] G'''' {r, r' , E) = 5^{r-r') . 

The numerical method for solving this equation is pre- 
sented in [l3|. At large energies, away from the threshold, 
the non-relativistic treatment is clearly not applicable. 
The spectral function in this domain, however, can be 
calculated using perturbation theory |l3| . 

To calculate G"''^ {r, r' , E) we need to specify the po- 
tential V{r) in the Schrodinger equation. A Coulomb 
plus linear form, phenomenologically successful in de- 
scribing quarkonia spectrum at T = 0, also gives a good 
parametrization of the lattice data on the static quark- 
antiquark potential for 0.1 fm < r < 1 fm. Only at 
distances r < 0.1 fm, not relevant for quarkonia stud- 
ies, the effect of the running coupling is important [l3{ . 
At large distances, rmed — 1-1 fm, the linear growth of 
the potential stops due to string breaking lj|. At high 
temperatures we expect that the effective range of the 
potential will be reduced and the interactions are expo- 
nentially screened at large distances. At sufficiently short 
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distances, on the other hand, the interaction between 
the heavy quark and antiquark is temperature indepen- 
dent Motivated by the above, we use the following 
parametrization for the potential 



or, 



V{r,T) 



Voo{T) 



a'jT) -u(T)r 



ed{T) , 
ed{T) . 
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Following the recent 2+1 flavor lattice QCD analysis [If 
we fix a = 0.385 and a = 1.263/rg , with ro = 0.469 fm 
being the Sommcr-scale determined in [l6l | . At zero tem- 
perature we choose Vmed = 1-1 fm and screening mass of 
/i = 0.4 GeV. The later is motivated by the fact that the 
heavy-light meson pair is in the isosinglet channel. This 
choice then corresponds to 14o — 1-2 GeV. The potential 
is shown in Fig. [1] For calculations of the charmonium 
spectral functions we have included in the potential also a 
relativistic spin-independent correction, estimated to be 
—0.8a/ (m^r) [17]. The zero temperature potential with 
the above choice of the parameters gives a fairly good 
description of the quarkonium spectrum. 

At finite temperature we fix the parameters in ^ uti- 
lizing the information available from lattice QCD on the 
free energy of static heavy quark- antiquark pair. Free en- 
ergy calculations are done in pure gluodynamics, 3-flavor 
and 2-flavor QCD and preliminary results are 

also available in the physically relevant case of one heavy 
strange quark and two light quarks [l5l| (quark masses 
correspond to pion mass of about 200 MeV). These calcu- 
lations show that the free energy is temperature indepen- 
dent for distances r < 0.4 im/{T/Tc) , while for distances 
rT > 0.8 i t is exponentially screened, with a screening 
mass estimated to be /i(T) = 1.4(1) ^/TTNjj6g{T)T 
[l^ [igj . Here g{T) is the 2-loop MS gauge couphng 
at scale 27rT. It was found also that at infinite sep- 
aration of the heavy quark and antiquark the free en- 
ergy approaches a constant Fao{T) . At high tempera- 
tures perturbation theory tells us that Foo(T) behaves 
like —o?J'^T ( 2 0| flattice calculations confirm this be- 
havior [ll, EEIII [ill)- For T > Tc we assume that 
FiyD{T) = a/T — bT, where a and b are temperature- 
independent constants. The first term can be viewed 
as the contribution from a non-perturbative dimension 
2 gluon-condensate 2l|, while the second term is the 



perturbative entropy contribution. This parametrization 
successfully describes the lattice data on Foo{T) in pure 
gauge theory [2lj , and we checked that it works very well 
for T > l.lTc also in 2-1-1 fiavor QCD. In our numerical 
analysis we use the value Tc — 0.192 GeV determined in 
(2^ . Although the 2+1 fiavor lattice data on the free en- 
ergy of static quark-antiquark used here were obtained on 
coarse lattices with temporal extent Nt = 4, discretiza- 
tion errors in this quantity are small [isj . 

It is important to note, that since it contains an en- 
tropy contribution, the free energy itself is not the po- 
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FIG. 1: The two choices of the potential (denoted as set 1 and 
set 2) used in our analysis for 1.2Tc and l.GTc. We also show 
the free energy of static quark anti-quark pair calculated in 
2+1 flavor QCD on the lattice [Tst - 



tential. It can, nevertheless, provide some constraints on 
the parameters of In particular, the discussion above 
implies 0.4 fm/(T/Tc) < r^ed < 0.8/T , and due to the 
negative entropy contribution Voo{T) > Foo{T). We use 
Voo = a/T and ^ = 1.4^/l + Nf /6g(T)T . The values 
of a.'{T) are determined in [l§|. The above described 
choice for Vqo, a' and /i, together with a requirement for 
the smoothness of the potential, predetermines the value 
of Tjned- This turns out to be close to 0.4 fm/(T/Tc), the 
value determined in lattice simulations. This physically 
motivated potential, labeled as set 1, is shown for differ- 
ent temperatures in Fig [TJ For comparison we display 
also the lattice data on the free energies. Denoted as 
set 2 in Fig [1] an alternative choice, an upper bound on 
the finite temperature potential is shown: When fixing 
fmed — 0.8/T the smoothness of the potential determines 
the value of I4o (T) . It turns out that this value is close to 
Uoo{T), the internal energy of the isolated static quark. 
We know that Uoo{T) provides an upper bound on the 
possible value of 14© (T) , because it contains all possible 
interactions of heavy quark with the medium. The values 
of ^J-{T) and a'{T) were fixed as in set 1. The numerical 
results on the spectral functions are mostly determined 
by the value of Voo(T) and are insensitive to other de- 
tails of the potential (see discussion in [13] also on the 
stability of the numerical solution). 

In what follows we present numerical results for these 
two choices of the potential, and derive an upper bound 
for the dissociation temperatures. In Fig. [2] we show 
the spectral functions above the deconfinement temper- 
ature. The S-wave charmonium spectral function at 
T = l.lTc exhibits a resonance peak with very small 
binding energy of about 0.014 GeV. Here, and in what 
follows, we define the binding energy as the distance be- 
tween the peak position and the continuum threshold, 
Ehin = 2mc.h + Voo — M , M being the resonance mass. 
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Above l.lTc the charmonium spectral functions show no 
resonance-hke structures, meaning that all charmonium 
states are dissolved. Even though there are no reso- 
nances, the spectral function is strongly enhanced com- 
pared to the non-interacting case. This is also illustrated 
in Fig. [51 We would like to point out, that this thresh- 
old enhancement compensates for the dissociation of the 
states, and thus dramatic changes seen in the spectral 
function are not reflected in the correlation function (for 



a detailed discussion see [121 )■ Also, the strong enhance- 



ment in the threshold region is an indication that the 
heavy quark and antiquark remain correlated. In the case 
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FIG. 2: S-wave charmonium (upper panel) and bottomonium 
(lower panel) spectral functions at different temperatures. 

of bottomonium we see only the ground state above de- 
confinement, all other states are dissolved. Furthermore, 
as Fig. [5] shows, there is no significant change in the peak 
position of the ground state up to 2Tc. Even though 
seemingly the resonance structure persists to tempera- 
tures even higher than this, the binding energy of the 
state is significantly reduced. 

When the binding energy of a resonance drops below 
the temperature the state is weakly bound, and ther- 
mal fluctuations can destroy it by transfering energy and 
exciting the quark anti-quark pair into the continuum. 
The rate of this excitation, or equivalently the width 
of the quarkonium states, is determined by the bind- 
ing energy [23| . Therefore, in order to provide an upper 
bound on the dissociation temperature we need to esti- 
mate an upper bound for the binding energy. To do this, 
we calculate quarkonium spectral functions for the set 2 
potential, providing the maximum possible binding still 
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FIG. 3: An upper limit of the binding energy of different 
quarkonium states (top) and the quarkonium width (bottom). 
The open squares show the width of the IS bottomonium 
state multiplied by six for better visibility, which has been 
calculated in the limit of small binding. 



consistent with the lattice results on screening. We find 
that with this choice of the potential the S-wave char- 
monium spectral function has resonance-like structures 
up to ~ l.BTc. Furthermore, we also see resonance like 
structures in the bottomonium spectral functions corre- 
sponding to the IP and 25* states. In the upper panel 
of Fig. [3] we show the corresponding binding energies of 
the different quarkonium states. Let us note, that in the 
past quarkonium widths at finite temperature have been 
calculated using perturbative QCD and the Boltzmann- 
approximation, assuming an ideal quark-gluon plasma. 
See '24] for a recent analysis. For quarkonium sizes re- 
alized in nature the validity of the perturbative calcula- 
tions of the quarkonium-gluon cross section is doubtful. 
Furthermore, the Boltzmann-approximation breaks down 
if the binding energy is smaller than the temperature. 
In ^231 the quarkonium dissociation rate due to thermal 
activation into the continuum has been estimated non- 
perturbatively, using a resonance plus a continuum model 
for the spectral function. The thermal dissociation rate 
r(T) has a particularly simple form in two limits [l^]: 
The limit of large and small binding, respectively: 

Y{T) = (^Afe-^--/^ , Eun » T , (4) 

r(T)-fvS' ^^™«^- (5) 

Here L is the size of the spatial region of the potential. 
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given by the distance from the average quarkonium ra- 
dius to the top of the potential, i.e. L = Vmed — (r^)^/^. 
From the top panel of Fig. [3]it is clear that for T > l.lTc 
all quarkonium states have binding energy smaller than 
the temperature, with the exception of the T(IS') state. 
Their width can thus be estimated using ([S]). The results 
are shown in the lower panel of Fig. [3l Note that all of 
these states have widths larger than 200 MeV, and are 
therefore likely to be dissociated in the plasma. In the 
case of the T(15') for T < l.QTc we use dU to estimate 
the width, which is found to be smaller than 40MeV. 
This is in fairly good agreement with the perturbative 
estimate of [2J]. For temperatures T > 1.6Tc even the 
IS" bottomonium is in the regime of small binding, and its 
width becomes very large by 2.6Tc (see Fig. [3|). Note that 
uncertainty in the value of Uoo in the lattice calculations, 
indicated as a band in Fig. [1] leads to uncertainty in the 
binding energy estimate of about 10%. When the ther- 
mal width is significantly larger than the binding energy 
the resonance structure seen in our calculation will not 
be observable in reality. We define the dissociation tem- 
perature as the smallest temperature where no resonance 
structure can be seen in the spectral function. The upper 
limit for the dissociation temperatures of the quarkonium 
states we determine by posing the conservative quantita- 
tive condition T{T) > 2E\jin{T). The corresponding dis- 
sociation temperatures are summarized in Table . A less 
conservative criterion r(r) > Ei,in{T) would reduce the 
dissociation temperature by roughly 10%. 



state 




^' 






Xb 


T 




< Tc 


< Tc 


1.2Tc 


1.2Tc 


l.3Tc 


2Tc 



TABLE I: Upper bound on dissociation temperatures. 



In conclusion, we determined quarkonia spectral func- 
tions in the quark-gluon plasma using a potential model 
with two choices for the potential, both motivated by lat- 
tice QCD results on the free energy of a static quark anti- 
quark pair. We found that, due to color screening, for 
the first chosen potential most quarkonia states, except 
the T, dissolve at temperatures close to that of decon- 
finement. For the most extreme potential which is still 
compatible with lattice data, resonance structure in the 
spectral functions exists up to higher temperatures. This 
potential provides an upper limit on the binding energy. 
Using the binding energy we calculate the width of var- 
ious states, and give upper bounds on their dissociation 
temperatures which are significantly lower than previous 
estimates. As such, the model proposed in [25[, where 
J/ip suppression is due only to melting of the Xc and 
rp' states, cannot explain the nuclear modification factor 
Raa measured in the experiments since color screening 
dissolves the J/tp- On the other hand, the enhancement 
of the spectral function near the threshold shows that 



the heavy quarks and antiquarks remain strongly cor- 
related in the plasma even though they do not form a 
bound state. This correlation could lead to the regenera- 
tion of some quarkonium states when the plasma converts 
to hadronic matter increasing Raa values above expec- 
tations from screening alone. The quark and antiquark 
may even reform into a higher excited state. For a quan- 
titative description of Raa, a model calculation of regen- 
eration effects is needed. More precise calculations of the 
spectral function and detailed lattice calculations of the 
static quark-antiquark correlators will also be helpful. 
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